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ABSTRACT 

The extinction curve derived from ISO mid-infrared (IR) observations of the 
Galactic center (GC) exhibits a surprisingly flat behavior in the ~ 3 to 8 /im 
region, contrary to the deep minimum expected from standard interstellar dust 
models consisting of bare silicate and graphite dust particles. We show that this 
extinction is likely caused by the presence of metallic needles in the interstellar 
medium (ISM) towards the Galactic center. If the needles contribute only to the 
3-8 /im extinction, they must have a long wavelength cutoff at ~ 8 /im, and 
therefore a typical length over radius ratio of ~ 600, smaller than the ~ 3 x 10 3 
aspect ratio determined for the needles in Cas A. Homogeneously distributed 
throughout the ISM, they comprise only a minor mass fraction of the ISM, with 
a needle-to-H mass ratio of ~ 5x 10~ 6 , which is equivalent to 0.14% of the 
silicate dust mass. Their total ISM abundance can then be readily explained 
by the combined production in SNe and O-rich stellar outflows. It is currently 
unclear how universal the GC extinction law is. Local 2-5 extinction 
measurements seem to be consistent with the standard extinction law, suggesting 
a non uniform distribution of needles in the ISM. The GC observations show that 
metallic needles, in spite of their low abundance or non uniform distribution, can 
be the dominant source of opacity in the 3-8 /im wavelength region. However, 
expelled into the intergalactic medium, their abundance is too low to cause any 
dimming of cosmological sources, and their length is too short to make them a 
significant source of submillimeter emission. 
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1. INTRODUCTION 

Analysis of the ~ 2.6 to 9 /zm spectrum of the Galactic center (GC) obtained with the 
Short Wavelength Spectrometer (SWS) on board the Infrared Space Observatory (ISO) by 
Lutz et al. (1996) and Lutz (1999) has revealed an anomalous extinction law, compared 
with that expected from standard interstellar dust models consisting of a population of bare 
silicate and graphite grains and polycyclic aromatic hydrocarbons (PAHs) (Li & Draine 
2001; Zubko, Dwek, & Arendt 2004, hereafter ZDA). The GC extinction is characterized by 
a relatively flat behavior at ~ 3 - 9 /im, contrary to model predictions of a deep minimum 
in that wavelength region. The minimum results from the steep decline in the graphite and 
silicate extinction from UV-optical to near-IR wavelengths, which is followed by a sharp rise 
at ~ 10 fim caused by the silicate 9.7 /im absorption feature (see Figure 1). ZDA have shown 
that an interstellar dust model consisting of bare silicate and amorphous carbon dust, PAHs, 
and composite particles containing silicates, organic refractories, water ice, and voids, can 
reproduce the GC extinction curve (see model COMP-AC-x in Figure 21 of ZDA). Here we 
pursue an alternative explanation, namely that the "anomalous" extinction is produced by 
a population of metallic needles in the interstellar medium (ISM). 

Metallic needles (whiskers) play an important role in Hoyle-Narlikar cosmologies. Their 
presence in the intergalactic medium has been invoked by Hoyle & Wickramasinghe (1999) to 
create the cosmic microwave background (CMB) by the thermalization of starlight. Needles 
can, in principle, provide the needed opacity at high redshift (Wright 1982), however, Li 

(2003) recently raised questions about the applicability of the antenna theory to modeling 
the optical properties of exceedingly thin needles, and the amount of intergalactic iron needed 
to produce the CMB. Banerjee et al. (2000) suggested that intergalactic absorbing needles 
can explain to the observed redshift-magnitude relation of Type la supernovae (Riess et al. 
1998, Perlmutter et al. 1999) without the need to invoke a positive cosmological constant. 
More recently, Narlikar et al. (2003) suggested that needles can also explain the CMB 
anisotropy detected by the Wilkinson Microwave Anisotropy Probe (Bennett et al. 2002). 

Here we suggest a more modest role for needles as a minor dust constituent in the 
Galactic ISM that nevertheless may be the dominant source of extinction in the ~ 3 to 
8 jum wavelength region. Conductive metallic needles have recently been proposed by Dwek 

(2004) as the source of the submillimeter emission detected with the SCUBA (Dunne et 
al. 2003) from Cas A. This, and the possible detection of needles in the ejecta of SN1987A 
(Wickramasinghe & Wickramasinghe 1993) may suggest that needles can be readily created 
in supernovae. Furthermore, Kemper et al. (2002) have suggested that non-spherical metallic 
particles can provide the 3-8 /xm opacity needed to fit the spectral energy distribution of 
the OH/IR star OH 127.8+0.0, suggesting that metallic needles may also be produced in 
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Fig. 1. — Extinction predicted by the BARE-GR-S dust model of Zubko, Dwek, & Arendt 
(2004; ZDA) with the additional population of metallic needles is compared to the Galactic 
center extinction derived by Lutz (1999). The various contributors to the extinction include 
PAHs (dashed-dotted line), graphite (dotted line), silicates (short dashed line), and metallic 
needles (long dashed line). The standard extinction predicted by the BARE-GR-S model is 
shown as a thin solid line, and falls short of the observations. The thick solid line shows is 
the total extinction obtained with the addition of needles to the population of interstellar 
dust particles. 
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quiescent O-rich stellar outflows. Conductive needles may therefore be a noticeable dust 
constituent of the general interstellar medium (ISM). 

Needle properties required to explain the GC extinction curve are derived in §2, and the 
astrophysical implication, whether the GC extinction is a special case or more characteristic 
of the general ISM is discussed in §3. 



2. NEEDLE PROPERTIES 

The extinction optical depth r(A) at wavelength A to a source located at distance S due 
to an assembly of dust particles along the line of sight (LOS) is given by: 
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where the sum is over dust composition j. The integral in square brackets is over the grain 
size distribution fj(a) (which is normalized to unity over the {a m i n j, a max ,j} size interval), 
rrid t j(a, s) = Anpjd 3 /3 is the mass of a dust particle of radius a and mass density pj located 
at s, and /t(A, a) is the mass extinction coefficient of the dust. The outer integral is over 
the LOS, where nn(s) is the H-number density at distance s. 

Equation (1) can be written in a simplified form as: 

r(A) = J2^)M d ALOS) (2) 

j 

j 

where Kj(\) represents an average over the grain size distribution, and M.d,j is the mass 
column density of the dust of composition j, pdj is the average dust mass density in the 
ISM, and L is the total distance along the LOS. 

Figure 1 depicts the average interstellar extinction curve for a mixture of bare graphite 
and silicate particles and PAHs that satisfies the UV-optical extinction curve, the diffuse 
IR emission, as well as the solar interstellar abundances constraints (model BARE-GR- 
S in Zubko, Dwek, & Arendt 2004, ZDA). The extinction was normalized to unity at K 
(A = 2.2 pm), for which r(K)=3.47/l. 086=3.20 (Rieke, Rieke, & Paul 1989, Lutz et al. 
1996). At ~ 10 /iin the extinction is dominated by the silicate absorption feature which has 
a mass absorption coefficient k, s u(9.7 pm) = 3xl0 3 cm 2 g -1 . Also shown in the figure are 
the new improved extinction measurements derived from the SWS recombination line data 
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presented by Lutz (1999). These measurements reinforce the previous results (Lutz et al. 
1996) that the ~ 3 - 8 fim wavelength extinction is essentially flat, with an average value 
of ~ 0.6 x t(K), and lacks the deep ~ 7 /im minimum characterizing the standard bare 
graphite-silicate interstellar dust model. 

Metallic needles with appropriate physical parameters can fill in the minimum in the 
"standard" interstellar extinction curve and produce the observed extinction towards the 
GC. At wavelengths below the long wavelength cutoff Ao (see below) one can use Mie cal- 
culations for infinitely long cylinders to calculate their optical properties. Their extinction 
is dominated by absorption, and for thin needles with a resistivity pn = 1 x 10~ 18 , and 
length £ » a, where the radius a = 0.01 /im, the Mie calculations give a mass absorption 
coefficient that slowly rises from values of ~ 2 x 10 5 , ~ 5x 10 5 , and ~9x 10 5 cm 2 g _1 at 2, 5 
and 10 /im, respectively, to a value of ~ 3 x 10 6 cm 2 g" 1 at 30 /im, (Hoyle & Wickramasinghe 
1999). 

The GC extinction measurement provide important constraints of the needle properties. 
The Cas A needle were required to be efficient submillimeter emitters, with a long wavelength 
cutoff in their mass absorption coefficient at 400 /im (Dwek 2004). Figure 1 shows that the 
derived A > 8.0 /im extinction is in fairly good agreement with that predicted by the standard 
bare silicate-graphite model. If the near-IR extinction were to be fitted with Cas A type 
needles, then their rising extinction efficiency at longer wavelengths would produce a ~ 
20 /im extinction that exceeds the observed value by a factor of three. Since the needles 
cannot contribute significantly to the opacity at mid-IR wavelengths, the cutoff wavelength 
Ao should be less than about 8 /im. In fact, an excellent fit to the GC extinction is obtained 
with A = 8 /im, after which the absorptivity of needles drops off as A~ 2 . The cutoff 
wavelength, A , is given by: 

x 1 (l/ Q ) 2 ^ 

Ao = 2 pRC WF) ' (3) 

giving a value of tja = 580 for a typical needle resistivity of 1 x 10 -6 Q cm « 1 x 10~ 18 s. This 
value is smaller than the value of £/a = 3000 found for the needles in Cas A, suggesting that 
grain-grain collisions or other break up mechanisms play an important role in the evolution 
of needles in the ISM. 

Figure 1 shows the contribution of the various dust components to the Galactic center 
extinction. The total extinction without the needles (thin solid line) exhibits a minimum 
at ~ 7 /im. This minimum is filled by the opacity of the needles (long dashed curve). The 
total extinction which includes the population of needles is shown as a thick solid line in the 
figure. 
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3. DISCUSSION 

An important quantity is the abundance of the needles needed to give rise to the observed 
GC extinction. Given a lifetime for their destruction in the general ISM, their abundance 
can be used to estimate their required production rate in supernovae and in stellar outflows. 
Furthermore, if intergalactic needles originate from galactic outflows, their abundance in 
the ISM of galaxies can be used to set an upper limit on their abundance in the general 
intergalactic medium. This upper limit can, in turn, be used to set an upper limit on the 
intergalactic opacity that can be attributed to metallic needles. 

Using eq. (2) we can write the needles-to-silicate mass ratio as: 

Pndl 
Psil 

At A = 9.7 pm, K si i = 3xl0 3 cm 2 g -1 , and = 5.0xl0 5 cm 2 g _1 . Figure 1 shows that for 
needles r n di(9.7 jum) = 0.22x3.20 = 0.70, whereas for the silicates t s u(9.7 /im) = 3.1. 

If we assume that the needles are uniformly mixed with the silicate dust, so that L ni ji = 
L si i, we get that needles contribute only a minor fraction to the total mass of ISM dust 
which, normalized to the abundance of silicates, is given by: 

= 1.4 x 10~ 3 for needles uniformly mixed in the ISM (5) 

Psil 

The silicate-to-H mass ratio is 0.004 (model BARE-GR-S of ZDA), giving a needle-to-H mass 
ratio of 5.6xl0~ 6 . 

Alternatively, we can assume that the abundance of needles is comparable to that of 
the silicates, i.e. p n di ~ p s a, in which case equation (4) implies that the needles are inho- 
mogeneously distributed in the general ISM, and occupy a small fraction of the LOS to the 
Galactic center, given by: 

— — ~ 1.4 x 10~ 3 for a patchy distribution of needles with p n< u ~ p S ii (6) 

L s il 

The universality of the Lutz extinction law towards the Galactic center is still an unre- 
solved issue. Analysis of stellar data (Whittet 1988) show that the 1.6 to 5 pm extinction 
towards stars within a distance of 3 kpc is fairly uniform, and closely represented by a A~ L7 
extinction law, similar to the standard extinction curve in Fig. 1. A similar extinction law 
was derived by Rieke, Rieke, & Paul (1989) in the same wavelength region from observations 
of select GC sources. Extended to 8 pm, these observations will argue against the universal- 
ity of the GC extinction law. However, a major uncertainty in these extinction measurements 
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is the intrinsic spectrum of the background stars, specifically the amount of excess photo- 
spheric IR emission, which may be variable. Near-IR (3-8 pm) extinction derived from the 
measurements of line ratios may therefore be more reliable than those derived from stellar 
continua, an issue that can only be resolved by more observations of background sources with 
well understood stellar spectrum. After the initial submission of this paper for publication, 
Indebetouw et al. (2004) presented the 1.25 to 8.0 pm extinction derived from data obtained 
by the Spitzer along the two £ = 42° and 284° lines of sight in the Galactic plane. The two 
extinction laws were found to be consistent with the Lutz (1999) measurements towards the 
GC, strengthening the case for the universality of this law in our Galaxy. 

Whether the needles are uniformly or inhomogeneously distributed in the ISM, the upper 
limit to the total mass of needles in the Galaxy is about 3xl0 4 M Q , assuming an ISM mass 
of 5xl0 9 M (Sodroski et al. 1997). The average lifetime of interstellar dust is about (4- 
6) x 10 8 yr (Dwek 1998, Jones 2001), requiring a needle production rate of ~6x 10~ 4 M yr _1 . 
The SCUBA observations of Cas A suggest that about 10~ 3 M of metallic needles formed 
in the SN ejecta. If the Cas A yield is typical, and if a comparable production rate takes 
also place in quiescent O-rich stellar outflows, then a Galactic SN rate of about 0.03 yr -1 
will be sufficient to account for the mass of ISM needles. 

It is interesting to examine whether the needles can provide any significant extinction 
to cosmological sources. The optical depth to redshift z is given by: 



where a(u) is the needles' cross section [assumed to be flat in the z/ to u /(l + z) frequency 
range], n n (z) oc (1 + z) 3 is their number density in the intergalactic medium, K n is their mass 
absorption coefficient, p n is their present mass density, and Rh = c/H is the Hubble radius. 
The function E(z) = [(1 + z) 2 (VL m z + 1) — z(2 + z)^] 1 / 2 , where Vt m and Vt^ are, respectively, 
the matter density and the the cosmological constant, normalized to their critical values. 
Using standard ACDM cosmology with fl m = 0.27 and = 0.73, the optical depth to a 
hypothetical source at redshift z = 2 is given by r = 4.5 K n p n Rn- The dust-to-gas mass ratio 
of needles in the Milky Way ISM sets an upper limit on their abundance in the intergalactic 
medium. The baryonic mass density for H = 70 km s _1 Mpc -1 is 4.2xl0~ 31 g cm" 3 , giving 
an intergalactic needle mass density of p n = 2.9xl0~ 36 g cur 3 . The Hubble radius Rh = 
4290 Mpc, giving a mass column density of ~ 2xl0~ 7 g cm" 2 out to z — 2. For a value 
of «(V) Rj 10 5 cm 2 g" 1 , the visual optical depth to that distance will be about 0.02. The 
needles giving rise to the GC extinction have therefore a negligible effect on the dimming 
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of light from cosmological sources, in particular high- ,2 Type la supernovae. They are also 
relatively short compared to the needles in the Cas A ejecta, and will therefore not give rise 
to any significant submillimeter emission. 

In summary the anomalously flat 3-8 /xm extinction observed by Lutz et al. (1996) 
and Lutz (1999) towards the Galactic center and by Indebetouw et al. (2004) towards the 
t = 42 ° and 284° lines of sight in the Galactic plane can be attributed to the presence 
of relatively short metallic needles in the ISM. The universality of this extinction law, and 
consequently the distribution of these needles in the general ISM is still a subject of ongoing 
investigation, that could be resolved with further Spitzer observations. 
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